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Intrinsic point defects in aluminum antimonide 
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Calculations within density functional theory on the basis of the local density approximation are carried out to 
study the properties of intrinsic point defects in aluminum antimonide. Special care is taken to address finite-size 
effects, band gap error, and symmetry reduction in the defect structures. The correction of the band gap is based 
on a set of GW calculations. The most important defects are identified to be the aluminum interstitial Al*^, 
the antimony antisites Sb^i and Sb^, and the aluminum vacancy V Al ~ . The intrinsic defect and charge carrier 
concentrations in the impurity-free material are calculated by self-consistently solving the charge neutrality 
equation. The impurity-free material is found to be n-type conducting at finite temperatures. 

PACS numbers: 61.72.J-. 61.72.Bb, 61.82.Fk, 71.15.Mb 



I. INTRODUCTION 

Aluminum antimonide is receiving much renewed interest 
for applications ranging from ionizing radiation detection to 
microelectronics to optoelectronics. For gamma radiation de- 
tection, AlSb is particularly promising as a novel material en- 
abling high energy-resolution detection at room temperature 
due to its indirect band gap of 1 .6 e V, the high atomic number 
of Sb, and the potentially high electron and hole mobilities of 
up to several hundred cm 2 V _1 s _1 at room temperature^ 
Improved materials for high resolution, room temperature 
gamma radiation detectors are critically needed for applica- 
tions in nuclear nonproliferation and monitoring, homeland 
security, and also medical and space imaging applications. 
Such detectors operate by counting the number of electron- 
hole pairs created in the semiconductor upon interaction with 
a gamma ray. Thus, a small band gap is generally desired 
to maximize the number of generated carriers, increasing the 
signal and reducing the shot noise. However, if the band gap 
is too small, excess thermal noise is generated from carriers 
thermally excited across the gap. At room temperature, a gap 
of ~1.6 eV is nearly optimal. In a similar vein, low intrin- 
sic carrier concentrations in the material may be desired to 
reduce the background signal, another noise source. Further- 
more, high energy resolution (counting statistics) is achieved 
by maximizing the efficiency of charge collection, which re- 
quires high carrier mobilities and long carrier lifetimes. An 
indirect band gap can be advantageous for maximizing carrier 
lifetimes by quenching radiative recombination. Finally, high 
atomic number in the material is desired to increase the stop- 
ping power for high energy radiation, reducing the required 
size of the device. Presently, the purity of large single-crystal 
growths of AlSb limits its performance for radiation detection 
application. 

In microelectronics, AlSb and related alloys containing In, 
are finding use in advanced field-effect transistor designs that 
promise higher switching speeds and lower power consump- 
tion compared to silicon devices^ The material has also 
found use in high current density, high speed resonant tun- 
nel diodes^ In optoelectronics, thin films of AlSb are being 
used in active regions and superlattice structures for claddings 
in novel type-II infrared cascade lasers^— Such lasers, emit- 
ting wavelengths around 3^4.3 /im, find application, for ex- 



ample, in remote atmospheric chemical sensing. Aluminum 
antimonide is particularly interesting in these applications be- 
cause it has a large 2.1 eV conduction band offset with InAs, 
which often comprises the other key component in the active 
layers of these devices. 



In the present work, we present a careful analysis of the 
thermodynamic and electronic properties of intrinsic point de- 
fects in AlSb, since the electronic and transport properties of 
the material can be degraded by detrimental defects. The main 
goals of this study are to identify the most important intrinsic 
point defects and to establish the concentrations of defects and 
charge carriers in thermal equilibrium. This work is part of a 
larger effort to understand the fundamental microscopic lim- 
its of performance of this and other semiconductor materials. 
Future work will focus on extrinsic impurities in the material, 
as well as the implications on carrier transport properties. 



Our major findings in this work are that the dominant native 
defects in AlSb are aluminum interstitials, antimony antisites, 
and aluminum vacancies, dependent on chemical environment 
and doping. The equilibrium concentration of total native de- 
fects near the melting temperature is found to be in the 10 16 
to 10 17 cm~ 3 range, while at lower temperatures, concentra- 
tions down to 10 10 cm~ 3 or lower are expected. The electron 
chemical potential in pure material is near the middle of the 
band gap, however the material tends naturally to be slightly 
n-type doped by charged aluminum interstitials. 



The paper is organized as follows. In Sec. [TT] we review 
the thermodynamic formalism we use to derive defect forma- 
tion energies and concentrations from first principles calcula- 
tions, and describe the computational details which underlie 
the present work. In Sec. Q]]] we report the relaxed defect ge- 
ometries, defect formation energies, and the charge carrier and 
defect concentrations, which are obtained by self-consistently 
solving the charge neutrality condition. Finally, in Sec.lIVIwe 
discuss the intrinsic limitations of AlSb on the basis of our 
results and the relation to experimental data. The Appendix 
contains a brief derivation of the band gap correction scheme 
utilized in the present work. 
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II. METHODOLOGY 



B. Computational details 



A. Point defect thermodynamics 

A material in thermodynamic equilibrium must contain a 
certain number of point defects at finite temperature, due 
to entropy. The temperature dependence of the equilibrium 
defect concentration can be shown to obey the following 
relation^ 



Cq exp 



AG D 
k R T 



(1) 



where cq denotes the concentration of possible defect sites, 
ks is Boltzmann's constant, and T is absolute temperature. 
The Gibbs energy of defect formation, AGd, can be split into 
three distinct terms 



AG D = AE D + TAS D + P AV L 



(2) 



The formation entropy, ASd, is typically on the order of 1 ks, 
so the entropy term hardly exceeds 0. 1 eV even at elevated 
temperatures. In some cases, entropic effects could play a 
role in stabilizing a defect at high temperatures when the en- 
thalpy differences are small. We do not explicitly treat this 
effect here, but point out cases in our results where it may be 
important. The formation volume, AVd, describes the pres- 
sure dependence of the Gibbs energy of formation and is typ- 
ically some small fraction of the atomic volume, so it can be 
neglected at ordinary pressures. The most important contribu- 
tion is the defect formation energy, AEp. In the following, 
we therefore focus on the formation energy and for simplicity 
ignore both the entropy and volume terms. 

For a binary compound the formation energy for a defect in 
charge state q is given b y 10 ' 11 

AE D =E d -^ (n A i + n Sb ) MaisI, 

- \ (tiai - n sb ) ( M b A f - /^f ) 

- \ ( n Ai - nsb) Afj, + q (£Vbm + Me) , (3) 

where Ejj is the total energy of the system containing the de- 
fect, n, denotes the number of atoms of type i, /j^ ulk is the 
chemical potential of component i in its reference state, and 
we have written the terms using explicit labels for AlSb. Ne- 
glecting entropic contributions, the chemical potentials of the 
reference phases can be replaced by their cohesive energies 
at OK. The formation energy depends on the chemical en- 
vironment via the parameter A/Lt, which describes the vari- 
ation of the chemical potentials under different conditions. 
The range of Afi is constrained by the formation energy of 
AlSb by \Afi\ < AH^ lsh , where for the present convention 
Afi = -AHf lsh and A^i = +AHf sh correspond to Al and 
Sb-rich conditions, respectively. Finally, the formation energy 
also depends on the electron chemical potential, fj, e , which is 
measured with respect to the valence band maximum, £Vbm- 



The energy terms in Eq. (0 were calculated using den- 
sity functional theory (DFT) carried out in the local den- 
sity approximation (LDA) using the Vienna Ab-initio Simu- 
lation Package (VASP)i2r— and the projector augmented-wave 
(PAW) method J^ii Defect formation energies were obtained 
using supercells of various size containing 32, 64, 128, and 
216 atoms. Extrapolation was used to account for finite-size 
effects as described in detail in Sec. Ill CI Brillouin zone in- 
tegrations were performed with fc-point grids generated using 
the Monkhorst-Pack schemed For the 32 and 64-atom cells, 
a non-shifted 6x6x6 mesh was used, while for the 128-atom 
cell, a shifted 3x3x3 grid was used. For the 216-atom cell, 
a non-shifted 4x4x4 mesh was constructed. The plane wave 
cutoff energy was set to 300 e V and Gaussian smearing with 
a width of 0. 1 eV was used to determine the occupation num- 
bers. For charged defect calculations, a homogeneous back- 
ground charge was employed (by omitting the G = term in 
the potential) to ensure charge neutrality of the entire cell. 

Atomic relaxations were performed to determine the equi- 
librium structures of the defects, with ionic forces converged 
to 20meV/A and all calculations performed at the theoretical 
equilibrium volume. Relaxations from various randomized 
initial configurations were performed to avoid high-symmetry 
local energy minima in the structures. 



C. Finite-size corrections 

In the supercell approximation there are spurious interac- 
tions between defects and their periodic images which lead to 
systematic errors^^^ For neutral defects the leading error is 
due to elastic interactions, which cause an overestimation of 
the formation energy. The strain energy of a point-like inclu- 
sion can be derived from linear elasticity theory and can be 
shown to fall off roughly with L~ 3 , where L is the distance 
between periodic imagesi2L2£ Therefore, the formation energy 
in the dilute limit (L — > oo) can be obtained by finite-size scal- 
ing with L~ 3 , which removes the elastic strain component. 

Makov and Payne considered the convergence of the en- 
ergy of charged systems in periodic systems and proposed a 
correction on the basis of a multipole expansion^ The lead- 
ing term corresponds to the monopole-monopole interaction 
and scales with L _1 . This term can be analytically deter- 
mined if the static dielectric constant of the medium, e, and 
the Madelung constant of the Bravais lattice of the supercell, 
a, are known;22i2i 



AE„ 



q 2 a 
~2Le' 



(4) 



The next higher order term in the expansion is the monopole- 
quadrupole interaction which scales as L~ 3 . Even higher 
order terms (0(L~ n ), n > 5) are usually small and there- 
fore neglected. In the present work, we have applied the 
monopole-monopole correction term using the experimental 
value for the static dielectric constant (e = 12). Then, since 
both the elastic and the monopole-quadrupole interactions 



3 



Number of atoms in ideal cell 
128 64 32 



3.5 



1.0 



R — n n _ 


- 






v*- 

Al 


B 




■ 


- «_ Al ' 


no mp-mp corr. 










~ ■ ^ 






-&- 


— e — 


© — — ^ — ~ _ ^ — 






Al 




'.Al 










—A 


A 








- A- 







0.0 0.5 1.0 1.5 2.0 2.5 3.0 3.5 
Concentration (%) 

FIG. 1. (Color online) Illustration of finite-size scaling of forma- 
tion energies to infinite dilution (concentration — > 0%) for three 
of the most important defects, with no band gap corrections and 
Afj, e = eV. The small filled and large open symbols respectively 
show the data without and with the monopole-monopole correction 
term. Note that, particularly for the charged aluminum vacancy, the 
monopole-monopole correction is the dominant size-dependent term. 
The monopole-monopole correction does not apply for neutral de- 
fects. 

scale with L~ 3 , we employed finite-size scaling with L~ 3 to 
correct for these terms. This extrapolation scheme gave very 
small extrapolation errors as shown in Tables [TT] and [Ell] The 
results of the finite-size scaling procedure are illustrated in 
Fig. Q] for the most important defects. Figure Q] clearly illus- 
trates the effect of the monopole-monopole correction term 
which tremendously reduces the variation between the super- 
cells. The remaining higher-order variations are very well 
captured by the L~ 3 finite-size scaling. In addition, the po- 
tential alignment correction described in Ref. 19 is implicitly 
taken into account by our extrapolation schemed 

D. Band gap corrections 

The underestimation of the band gap (by the LDA) affects 
the formation energies, as discussed in detail in the Appendix. 
A simple correction scheme based on the band energy has 
been proposed by Persson et al^> and is further motivated in 
the Appendix. The correction amounts to the following term 
which is added to the as-calculated formation energies: 

AE corr = ( q+ Azh ) AEvB + Az e AE CBl (5) 

where AEVb and AEcb are shifts of the valence and conduc- 
tion band edges, respectively, to correct the band gap, while 
Az/j and Az e are the number of unoccupied valence and oc- 
cupied conduction band states, respectively. 

In the present work the energy offsets AEVb and AEcb 
were obtained from GqWq calculations— within the single 



plasmon-pole model as implemented in ABlNlTias— Non- 
selfconsistent Go Wo calculations were employed to properly 
refer the quasiparticle energies to the same potential zero as 
the LDA eigenvalues. Fritz-Haber-Institute norm-conserving 
pseudopotentials 31 in the Troullier-Martins scheme^ were 
used with a cutoff energy of 15 Hartree (Ha). The other rele- 
vant cutoff energies used in the calculation were 5 Ha for the 
self-energy wave functions, 6 Ha for the exchange part of the 
self-energy, and 6 Ha for the screening matrix. The number 
of bands in the self-energy and screening matrix calculations 
were 100 and 150, respectively. 

In Tablelnlwe report both the as-calculated formation ener- 
gies (including finite-size scaling) and the band gap-corrected 
formation energies. The correction terms can be reproduced 
using the values for Az/, and Az e included in the table. Note 
that we extracted Azh and Az e values only for the most im- 
portant defects for which an unambiguous distinction between 
valence and conduction band states was possible. A direct 
comparison of the as-calculated and band gap-corrected val- 
ues is shown in Fig. [3] and discussed later in Sec. IIII Cl 

The defect concentrations calculated in Sec. IIII Dl were ob- 
tained using the band gap-corrected formation energies and 
the GW band gap. 



III. RESULTS 
A. Bulk properties 

As described in Sec. Ill Al and evident in Eq. (O, the deter- 
mination of defect formation energies requires that we also 
calculate bulk properties of the solid and its constituents in 
reference states. The reference states for Al and Sb are face- 
centered cubic (fee) and rhombohedral solids, respectively. 

For fee aluminum, we obtain a lattice constant of 3.99 A 
and a cohesive energy of —4.19 eV/atom which compare rea- 
sonably well with the experimental values of 4.05 A (room 
temperature )22 and —3.38 eV/atom^i The underestimation of 
the lattice constant and the overestimation of the cohesive en- 
ergy are typical for LDA calculations. Antimony has a rhom- 
bohedral groundstate structure (R3mh, space group no. 166, 
Strukturbericht symbol A7) for which the DFT calculations 
yield a lattice constant of 4.46 A and a rhombohedral angle 
of 59.0° (experimental values: 4.50 A and 57.1°, Ref. Q) 
and a cohesive energy of —4.81 eV/atom (experimental value: 
-2.72eV/atom, Ref. (H). 

At ambient conditions bulk aluminum antimonide adopts 
the zinc -blende structure (F43m, space group no. 216, Struk- 
turbericht symbol B3). The DFT calculated lattice constant 
is 6. 12 A in good agreement with the experimental value 
of 6.13 A (300 K). The DFT calculations furthermore yield 
a formation enthalpy of — 0.28eV/f.u. (experimental value: 
— 0.84eV/f.u.). The direct and indirect band gaps are calcu- 
lated as 1.53 and 1.12eV, respectively, while the experimental 
values are 2.30 and 1.62 eV at 300 

From the GW calculations, we obtain direct and indirect 
band gaps of 2.3 1 and 1 .64 eV, respectively, in excellent agree- 
ment with the experimental values cited above. The shifts of 
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FIG. 2. (Color online) Relaxed defect geometries for (a,b) vacancies, (c,d) antisites, (f,g,h) several tetrahedral interstitials, and (i) one split 
interstitial, (e) The conventional 8-atom unit cell is shown for reference. Purple (dark grey) and beige (light grey) balls represent antimony and 
aluminum atoms, respectively. Note that the negatively charged antimony vacancy in (b) shows a pronounced Jahn-Teller distortion. (Indicated 
bond lengths were taken from 216-atom supercells.) 



TABLE I. Stability of split interstitial configurations. 
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' Relaxations starting from the idealized positions maintained the initial symmetry, but when started from 
randomized positions, relaxed to the indicated structure 
* Split interstitial displaced along (100) from ideal position 



the valence and conduction band edges showed small varia- 
tions with fc-vector, so we calculated Ai?vB and AEcb which 
appear in Eq. (0 as weighted averages over all fc-points 
included in the calculations. The shifts thus obtained are 
AE YB = -0.44 eV and AEcb = 0.16eV. 



B. Defect structures 

We considered all possible native defects in AlSb up to split 
interstitials. In total, this amounts to 18 different defects (two 
vacancies, two antisite defects, four tetrahedral interstitials, 
two hexagonal interstitials, and eight split interstitials), not 
all of which are stable. For each defect, we investigated a 
series of charge states, generally from q = —3 to q = +3, as 
appropriate. Defect complexes were not considered. 

Vacancies. In the ideal zinc-blende structure, both Al and 
Sb sites possess tetrahedral symmetry with nearest neighbor 
distance of 2.65 A. If one removes a single atom and allows 
the system to relax from randomized positions, the aluminum 



vacancy Vm maintains the Td symmetry for all relevant charge 
states and the surrounding antimony ions relax inward by 
0.36 A (for q = 0) to 0.38 A (for q = -3) [see for exam- 
ple Fig. |2(a)| . This behavior is typical of cation vacancies in 
III-V and II- VI zinc -blende semiconductors, where a dimer- 
ization transformation of the atoms surrounding the vacancy 
is generally not energetically favorable^ In contrast, the an- 
timony vacancy Vsb exhibits a Jahn-Teller distortion to a local 
tetragonal symmetry for all but the positive charge states. The 
relaxed configuration for is shown in Fig. |2(b)| represen- 
tatively, also indicating the pairing of Al atoms in the first 
neighbor shell of the vacancy. 

Antisites. Both the A1sh and SbAi antisite defects maintain 
the Td symmetry for all relevant charge states. In the case 
of Alsb, the surrounding aluminum ions relax inwards, while 
for SbAi the nearest neighbor antimony ions relax outwards 
[compare Figures |2|c) and (d)]. These relaxations occur as 
expected based on the atomic radii. 

Tetrahedral interstitials. In the zinc -blende structure there 
are two distinct tetrahedral sites: one centered on an Al tetra- 
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FIG. 3. (Color online) Dependence of the formation energies on the electron chemical potential under Al-rich (left, A/x = —AHf), stoi- 
chiometric (center, Ap = eV), and Sb-rich (right, A/j, = +AHf) conditions for the most important point defects. The top panels show the 
formation energies based on the as-calculated values, while the bottom panels show the results of applying the band gap correction term of 



Eq. {3). The slope of each line is determined by the charge state, according to Eq. 10. (AH'? < 



-0.28 eV) 



hedron (4d site) and one centered on an Sb tetrahedron (4b 
site). Thus, there are four possible types of tetrahedral inter- 
stitials: AI^.ai and Sb^At on the 4d site; Al^sb and Sb,,sb on 
the 4b site. 

Both kinds of tetrahedral aluminum interstitials (Al,,*!, 
Ali sb) maintain the Td symmetry after relaxation, with the 
neighboring ions relaxing outwards as shown in Figs.|2jf) and 
(g). The Sbi.sb configuration is unstable in all but the +2 
charge state; in all other charge states, the interstitial atom 
relaxes either onto a hexagonal site or forms a split intersti- 
tial. Conversely, Sb^Ai is stable in all of its charge states, but 
does exhibit Jahn-Teller distortions as illustrated in Fig. |2(h)| 

Hexagonal interstitials. The hexagonal interstitial is lo- 
cated at Wyckoff position 16e. The aluminum hexagonal in- 
terstitial Ali hex was found to be unstable for all charge states, 
with the Al atom observed to relax into the Ali Ai configura- 
tion even when starting from ideal positions. For the antimony 
hexagonal interstitial Sb^hexi the +2 charge state was found 
to relax directly into the Sb^Ai position, whereas for charge 
states +1 and —2 the interstitial atom remained in the hexag- 
onal site. The other charge states when perturbed from the 
ideal position relax into the (Sb — Sb)sb(no) split interstitial, 
described below. 

Split interstitials. Split interstitial configurations, in which 
two atoms share one atom site, have been extensively dis- 
cussed in the literature for the zinc-blende structurejii^^ 



We have carried out an exhaustive exploration of the struc- 
tures of these defects in AISb. The following cases were 
considered: two Al atoms oriented along (100) sharing one 
Al site [(Al — AI)ai(ioo)]' me same atoms but oriented along 
(110) [(Al — Al) A i(iio)L one Al atom and one Sb atom ori- 
ented along (100) sharing one Al site [(Sb — A1)ai(ioo)L 
and the same combination of atoms oriented along (110) 
[(Sb — AI)ai(ho)]- Four more configurations are obtained cor- 
responding to the same combinations above, but sharing an Sb 
site. 

Most of these configurations actually are found to be unsta- 
ble with respect to other interstitial configurations. The results 
of these calculations are summarized in Table U which shows 
which configurations were found to be stable and which ones 
were unstable or only conditionally stable with respect to al- 
ternative interstitial configurations. In certain cases, a struc- 
ture relaxed starting from the ideal atomic coordinates main- 
tained the starting symmetry after relaxation, but when started 
from randomized coordinates, the structure relaxed to a differ- 
ent configuration; these cases are indicated in Table Q]by list- 
ing both final configurations. The cases simply marked stable 
(single entry 'stable') in Table [Jrelaxed to the ideal symmetry 
configuration for all starting configurations. 

For completeness, we note that the split interstitial 
(Sb — Al) A i( 10 o) in charge states q = 0, —1, and —2 is dis- 
placed along (100) from the ideal position, such that the Sb 
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interstitial is located in between two regular Sb atoms along 
the (110) direction. 

C. Formation energies 

The defect formation energies calculated using Eq. © for 
Al-rich and Sb-rich conditions and an electron chemical po- 
tential at the valence band maximum are given in TablesHIland 
HTT1 The formation energies for the dominant defects (lowest 
AEjj) are shown as a function of the electron chemical po- 
tential in Fig. [3] 

Figure[3]also illustrates the effect of the band gap correction 
described in Sec. Ill Dl In the top row we show the formation 
energies after the finite-size scaling procedure was applied but 
without the band gap correction. The bottom row shows the 
results including the band gap corrections. It is apparent that 
as the band gap correction is applied donor and acceptor lev- 
els track the conduction and valence band edges, respectively. 
This feature is independent of the relative values of A£Vb and 
AE CB in Eq. ©. 

It is interesting to compare the band gap-corrected forma- 
tion energies in the bottom panel of Fig. [3] with those in the 
top panel if the band gap is simply extended to the experimen- 
tal value by shifting the conduction band upwards. Consider- 
able differences are observed between the two cases, with the 
values corrected using the scheme presented here being more 
consistent. In fact, we observe that, qualitatively, the results 
obtained by applying no correction at all are more similar to 
the corrected results than what is obtained by assigning the 
band gap error fully to the conduction band42iii Furthermore, 
a correction procedure often suggested in the literature^ to 
simply shift the ionization energies of donor-like defects to 
track the conduction band minimum (implicitly leaving the 
ionization energies of acceptor-like defects tracking the va- 
lence band maximum), without accounting for the occupation 
of states, does not properly correct the errors in the formation 
energies. The scheme employed here results in an identical 
shift of the ionization energies, but also corrects the errors 
in the formation energies. However, we note that this scheme 
still neglects both level relaxations and changes in the double 
counting term. 

We believe the method presented here is the most consis- 
tent way to address the LDA band gap problem, with the use 
of GW calculations providing a first principles approach to 
calculating the correction terms. We find that neglecting to in- 
clude the band gap correction terms in the formation energies 
leads to significant errors in the prediction of defect concen- 
trations and which defects are dominant. 



D. Defect and charge carrier concentrations 

With the formation energies known, the equilibrium defect 
concentrations for a given chemical potential difference A/i 
can be calculated using Eq. (|T). The defect concentrations 
depend on the electron chemical potential via Eq. ([3]). In the 
absence of extrinsic defects, the electron chemical potential is 



TABLE II. Formation energies, in eV, of intrinsic point defects un- 
der Al-rich (A^ = -AH f ) and Sb-rich (Aft = +AH f ) condi- 
tions for an electron chemical potential at the valence band maxi- 
mum [/i e = OeV in equation Eq. (3}]. Both the as-calculated and 
the band-gap corrected formation energies are given. The number 
of occupied (unoccupied) conduction (valence) band states defined 
through Eq. dA.4b is given in the third column (Az e ,j,), where neg- 
ative (positive) values correspond to electrons (holes). The extrap- 
olation error, in eV, of the finite-size scaling procedure is given in 
brackets in the last column. 
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2.66 


(0.01) 




+ 1 





1.15 


1.43 


0.71 


0.98 


(0.01) 




+2 


+ 1 


1.44 


1.72 


1.00 


1.27 


(0.02) 




"TO 




1 

1 .yo 
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z.zz 






^U.UD ) 


Ali.sb 


-2 


— 


5.49 


5.77 


— 


— 


(0.07) 




-1 


-2 


4.01 


4.28 


4.78 


5.06 


(0.03) 







-1 


2.75 


3.03 


2.91 


3.19 


(0.01) 




+ 1 





1.73 


2.01 


1.29 


1.56 


(0.01) 




-1-9 




1 .Do 


i.yo 


1 9/1 


1.01 


(n m ^ 




+3 




2.00 


2.28 






(0.01) 


VA3 


-3 





3.99 


3.72 


5.32 


5.04 


(0.02) 




-2 


+1 


3.35 


3.07 


4.68 


4.40 


(0.02) 




-1 


+2 


2.94 


2.66 


4.27 


3.99 


(0.03) 







+3 


2.75 


2.47 


4.08 


3.80 


(0.04) 


Sb A i 


-1 




3.25 


2.70 






(0.03) 










2.03 


1.48 


2.03 


1.48 


(0.02) 




+ 1 


+1 


1.52 


0.97 


1.52 


0.97 


(0.01) 




+2 


+2 


1.40 


0.85 


1.40 


0.85 


(0.01) 


Sbi,sb 


+2 




3.63 


3.35 






(0.01) 


Sbi.Ai 


-2 




6.32 


6.04 






(0.10) 




-1 


-4 


4.84 


4.57 


5.94 


5.66 


(0.02) 







-3 


3.92 


3.65 


4.41 


4.13 


(< 0.01) 




+ 1 


-2 


3.16 


2.88 


3.04 


2.77 


(0.13) 




+2 


-1 


2.93 


2.66 


2.21 


1.93 


(0.01) 


Sbjfoe^ 


-2 




6.24 


5.96 






(0.09) 




-1 


-4 


4.96 


4.68 


6.06 


5.78 


(0.02) 







-3 


3.93 


3.66 


4.42 


4.14 


(0.01) 




+ 1 


-2 


3.14 


2.86 


3.02 


2.75 


(0.01) 
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Sb— rich 



Al— rich 



TABLE III. As-calculated formation energies for the split-interstitial 
configurations which were found to be stable or metastable. Values 
are given for both Al-rich and Sb-rich conditions. The extrapolation 
error of the finite-size scaling procedure is given in brackets in the 
last column. Values are in eV. 



Defect 


q 


Al-rich 


Sb-rich 




(Al - A1)ai (10 o> 


-2 


6.93 


7.20 


(0.11) 




-l 


5.48 


5.76 


(0.08) 







4.27 


4.54 


(0.06) 




+ 1 


3.31 


3.58 


(0.06) 




+2 


3.38 


3.65 


(0.02) 


(Al - A1) A | (110> 


-2 


4.74 


5.02 


(0.04) 




-1 


3.32 


3.60 


(0.02) 







2.68 


2.96 


(< 0.01) 




I 1 

+ 1 


O Of? 

2.o6 


2.04 


(0.01) 




+2 


2.54 


2.82 


(0.01) 


(Al - Sb) sb <iio> 


-2 


5.68 


5.96 


(0.06) 




-1 


4.34 


4.61 


(0.04) 







3.34 


3.62 


(0.01) 




I 1 

+ 1 


2.08 


2.yo 


(0.01) 




+2 


1.68 


1.96 


(0.01) 


(Sb - Alj A1(100 ) 


— 2 


0.94 


0.0/ 


(0.11) 




-1 


4.53 


4.26 


(0.09) 




() 


3.91 


3.63 


(0.03) 


(Sb - A1) AI(110) 


-2 


5.70 


5.42 


(0.08) 




-1 


4.26 


3.98 


(0.05) 







3.46 


3.19 


(0.01) 




+ 1 


3.01 


2.74 


(0.07) 




+2 


3.00 


2.72 


(0.01) 


(Sb — Sb)sb(ioo) 


+ 1 


3.72 


3.44 


(0.02) 




+2 


3.72 


3.44 


(0.01) 


(Sb - Sb) sb(110) 


-2 


5.13 


4.85 


(0.11) 




-1 


3.72 


3.45 


(0.08) 







3.28 


3.00 


(0.05) 




+ 1 


3.17 


2.90 


(0.03) 




+2 


3.45 


3.18 


(0.04) 



constrained by the charge neutrality condition 

defects 

= n e - n h - 2J qiCi, 
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FIG. 4. (Color online) Concentrations of individual defects at a tem- 
perature of 1300 K as a function of the chemical potential difference 
Afi. Line thicknesses scale with charge state. Under Sb-rich condi- 
tions the dominant defect is the neutral Sb A! antisite, whereas under 
Al-rich conditions the positively charged aluminum interstitial Al* A1 
is dominant. The crossover point between Al] A1 and Sb A1 (open cir- 
cle) determines the minimum in the total defect concentration shown 
in Fig. |5(a)[ while the crossover point between Al* A1 and Sb^j 1 " (filled 
circle) determines the minimum in the net electron concentration 
shown in Fig. |5(b)| The horizontal bar marks the calculated range 
of variation of A/i given by the formation enthalpy of the compound 
(-AH f <A/j,< +AH f , with AHf' = -0.28 eV). 



The defect concentrations calculated using the band gap 
corrected formation energies are shown in Fig. [4] for a rep- 
resentative temperature of 1300K (T^ b = 1327 K), with the 
line thicknesses indicating the charge state \q\. Figure|5]shows 
the dependence of the total defect and net electron concentra- 
tions on the chemical environment (chemical potential differ- 
ence A/i) for a variety of temperatures. For all cases shown 
here, the intrinsic (self-consistent) electron chemical poten- 
tial is located near the middle of the gap, although there are 
slightly more electrons than holes in the material (intrinsically 
n-type material). 



(6) 



IV. DISCUSSION 



since the intrinsic concentrations of electrons and holes, n £ 
and rih, respectively, are given by 



n e = j D(E)f(E;f* e )dE 

n h = J D(E) [l-f(E;/i e )]dE. 



(7a) 



(7b) 



Here, D(E) is the electronic density of states and f(E; /i e ) = 
{l + e (- E -^e)/ fc s T } is the Fermi-Dirac distribution. The 
implicit dependence of the charge neutrality condition on the 
electron chemical potential /j, e is apparent from Eqs. ©. To 
obtain the charge carrier and defect concentrations, then, we 
must iteratively solve Eq. (O to self-consistently determine 
the intrinsic electron chemical potential. 



An inspection of the Go Wq -corrected results in the lower 
panels of Fig. [3]shows that four different defects have the low- 
est formation energy and are thus the most abundant, depend- 
ing on the values of the electron chemical potential, /i e , and 
the chemical potential difference, A/i. The aluminum tetra- 
hedral interstitial AvJ^ is the dominant defect for fi e in the 
lower half of the band gap (p-type material), while the alu- 
minum vacancy is dominant for fi e in the upper half of 
the band gap (n-type material). For Al-rich conditions (left 
panel in Fig. [3), the AlSb antisite defect can also be impor- 
tant when the material is lightly doped ?i-type (near the cross- 
ing with V^~), particularly considering the uncertainty in the 
range of A// alc vs Afi expt and the possible errors in the for- 
mation energies of up to about 0.1 eV from entropic effects. 
For Sb-rich conditions (right panel in Fig. [3j, the SbAi antisite 
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-0.4 -0.2 0.0 0.2 0.4 -0.4 -0.2 0.0 0.2 0.4 
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FIG. 5. (Color online) (a) Equilibrium total defect concentration and (b) net electron concentration (n e — nu) at different temperatures as 
a function of the chemical potential difference. Note that the minimum in the total defect concentration corresponds to the crossover point 
between Al^ and Sb A i, while the minimum in the net electron concentration corresponds to the crossover point between Al^ and Sb^, 
as indicated in Fig. [4] The dominant defects on the Sb-rich side are SbAi and on the Al-rich side are AI^ai- The horizontal bar marks the 
calculated range of variation of A/i given by the formation enthalpy of the compound (—AHf < A/i < +AHf, with AHf ]c = —0.28 eV). 



defect (neutral or positively charged) is important over a wide 
range of /j, e , from the p-type regime to near the middle of the 
gap. However, we note that for /i e near the middle of the gap 
(intrinsic or compensated material), the dominant charged na- 
tive defect is always Al^, since the relevant Sb^ defect is 
uncharged. 

In the absence of impurities, the intrinsic electron chemi- 
cal potential is located near the middle of the band gap (see 
Sec. IIII DY Under this condition, the two relevant defects are 
therefore the Al i ^ 1 interstitial and the SbAi antisite, as illus- 
trated clearly in Fig. |4] By comparison with Fig. |5(a)| we see 
that the transition point between Al.^, and Sb^i that occurs 
close to A/i = eV corresponds to the minimum in the total 
defect concentration. In Fig. |5(a)| the dominant defect on the 
Sb-rich (left) side is Sb^j, while on the Al-rich (right) side it 
is Al^j. Near the transition (minimum in the curve), the two 
concentrations are comparable. In contrast, the minimum in 
the net electron concentration in Fig. |5(b)| is located on the Sb- 
rich side, corresponding to the crossing of the concentrations 
of Al^! and Sb^+ in Fig. g] Since the charge neutrality con- 
dition [Eq. ©] depends only on the concentrations of charged 
defects, the neutral Sb^, antisites do not contribute to charge 
compensation or net electron concentration. 

As mentioned in Sec. IIII Dl the calculated electron concen- 
tration slightly exceeds the hole concentration for pure ma- 
terial, yielding n-type intrinsic material irrespective of the 
chemical environment characterized by AyU. This behavior re- 
sults because the formation energies of acceptor-type defects 
(Vj^~) always exceed the formation energies of donor-type de- 
fects (Al^j, Sb^) when the electron chemical potential is 
close to the middle of the band gap (see Fig. [3]). It should be 
noted that this situation is changed if the material is extrinsi- 
cally doped. For /9-doped material, the donor-type native de- 
fects remain dominant and partially compensate the extrinsic 
dopant. For n-doped material, the ordering of the formation 



energies is reversed, but the then-dominant acceptor-type na- 
tive defects still partially compensate the extrinsic dopant. 

The temperatures indicated in Figs. [4] and [5] refer to ther- 
mal equilibrium conditions at those temperatures. In practice, 
these temperatures can be interpreted as corresponding to an- 
nealing temperatures, with the chemical potential difference 
referring to the chemical environment of the annealing process 
(e.g., an Sb overpressure corresponds to Sb-rich conditions; 
conversely, growth is often performed under Al-rich condi- 
tions). The highest temperature considered here, 1300K, is 
just below the melting point of 1327 K and might represent 
melt growth conditions. However, growth is typically per- 
formed too rapidly to allow equilibrium to be achieved and 
the non-equilibrium grown-in defect concentrations will be 
higher than predicted here. The curves in Fig. [5] essentially 
represent the predicted concentrations for infinitely long an- 
neals at the specified temperatures. 

At 1300K, the calculated net electron concentration varies 
roughly between 10 16 and 10 17 cm -3 depending on the chem- 
ical potential difference, A/x. If all defects are assumed to 
be sufficiently mobile down to 700 K so that the material can 
reach thermal equilibrium at that temperature, then the lower- 
most curves in Fig. [5] predict a net electron concentration of 
10 10 to 10 11 cm~ 3 and a total defect concentration between 
10 10 and 10 12 cm -3 . Since the diffusivity depends exponen- 
tially on the inverse temperature, the defect mobilities de- 
crease sharply with temperature. Therefore, as the tempera- 
ture is further lowered the system will no longer be able to 
reach the equilibrium concentration in reasonable time, which 
requires excess defects either to diffuse to the surface or to an- 
neal by recombination. In contrast to lattice defects, the intrin- 
sic electron and hole concentrations, n e and rih, readily adjust 
to temperature changes. The "freezing in" of the defect con- 
centrations is therefore expected to crucially affect the charge 
neutrality condition upon cooling to low (e.g., room) tempera- 
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tures (in particular if there are no extrinsic dopants). However, 
since the diffusivities of the individual defects are currently 
unknown, we cannot quantitatively describe this "freezing in" 
of the defect distributions in the present study; thus, predic- 
tions of the charge carrier and defect concentrations near room 
temperature and below may be unreliable. The determination 
of diffusivities for specific defects, to account for the kinetics 
of defect "freeze-in," will be the subject of future work. 

Experimentally, AlSb crystals grown from the melt have of- 
ten been found to display /?-type conductivity. 45 The present 
finding that the pure material behaves intrinsically n-type is, 
however, not in contradiction with this observation. A typi- 
cal experimental setup employs a graphite susceptor, alumina 
crucible, and quartz tubes for melting Sb, which are poten- 
tial sources of various impurities, most importantly carbon, 
oxygen, silicon, and aluminum. 45 In particular, carbon impu- 
rities act as acceptors, therefore accidental p-type doping is 
a very likely scenario which is supported by chemical analy- 
sis of grown AlSb crystals and measurements on intentionally 
doped samples^ The measured levels of carbon impurities 
and the experimental observation that the material becomes 
intrinsic around 1000 K (Ref. l35lo are consistent with our cal- 
culated excess native electron concentration at that tempera- 
ture. A detailed study of the role of extrinsic defects in AlSb 
is beyond the scope of the present work; however, the results 
of ongoing work to elucidate the effects of impurities and to 
investigate ways to optimize the electronic properties of the 
material by intentional doping are forthcoming. 



V. CONCLUSIONS 

In summary, we have employed density functional theory 
calculations to study the properties of intrinsic point defects 
in aluminum antimonide. An exhaustive set of defect con- 
figurations — including vacancies, antisites, and interstitials 
(tetrahedral, hexagonal, and split), with all relevant charge 
states — was considered based on knowledge of other III-V 
compounds. Relaxed atomic structures of each defect were 
carefully determined, and formation energies were calculated 
to evaluate the equilibrium concentrations of each defect. 
Strain and electrostatic artifacts related to the use of the super- 
cell approach were carefully removed by employing a finite- 
size scaling procedure, which involved performing a series of 
calculations for each defect with supercell sizes ranging from 
32 to 216 atoms. The underestimation of the band gap due 
to the local density approximation was taken into account by 
applying an a posteriori correction scheme that utilized sep- 
arate GW calculations to obtain valence and conduction band 
offsets, which enter the correction. 

Aluminum interstitials (Al,^), antimony antisites (Sb^,, 
Sb^ 1 "), and aluminum vacancies (V^~) were found to be the 
most dominant defects, depending on the electron chemical 
potential and the chemical potential difference (chemical en- 
vironment) of the system. We observe that Al^ intersti- 
tials and SbAi antisites dominate under Al-rich and Sb-rich 
conditions, respectively. Calculated formation energies were 
employed in solving the charge neutrality condition to ob- 



tain self-consistent defect concentrations and intrinsic elec- 
tron chemical potential for the pure material at various tem- 
peratures. We find the material to be intrinsically weakly n- 
type and predict both the total defect and the net electron 
concentrations. Near the melting point, the equilibrium con- 
centration of native defects is predicted to be in the 10 16 to 
10 17 cm" 3 range, while at lower temperatures, it is expected 
that concentrations down to the 10 10 cm -3 range or lower can 
be achieved. The net excess electron density in bulk grown 
material might be as high as 10 9 to 10 11 cm -3 from "freeze- 
in" of defects from melt solidification. 

For extrinsically doped material, which we do not treat ex- 
plicitly in detail in this work, the dominant native defects 
depend on the nature of the doping. For n-doped material, 
and Alj|^ tend to be important, while for /j-doped ma- 
terial, Al^j and Sb^ 1 " are important, depending on chemi- 
cal environment (Al-rich vs. Sb-rich). Some amount of self- 
compensation from the native defects occurs in both cases. 

Finally, we note that the present work is part of a concerted 
research effort which ultimately aims to provide a complete 
and consistent picture of the point defect properties of AlSb 
and the relations to carrier transport properties. We are en- 
gaged in further theoretical and experimental work to explore 
the role of extrinsic defects, as well as the scattering behav- 
ior of defects on carrier transport. In this context, the present 
study forms the basis for these future studies, which will be 
the subjects of forthcoming reports. 
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Appendix: Band gap correction of defect formation energies 

While DFT calculations typically give very reasonable val- 
ues for energy differences within a group of bands (e.g., the 
valence band), energy differences between different groups of 
bands (i.e., band gaps, e.g., between the valence and conduc- 
tion bands) are much less reliable. This shortcoming particu- 
larly affects the differences between the valence and conduc- 
tion bands, giving rise to the well-known band gap error. 

The incorrect description of the energy differences between 
different groups of bands can affect the total energy. The most 
sensitive contribution is the band energy which is given by 

i k 
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where i and k run over bands and fc -points, respectively, and 
fik and Sik are the occupation numbers and eigenvalues, re- 
spectively. Without loss of generality, one can divide the band 
energy into separate sums over the valence and conduction 
band states, as 

VB CB 

E b = E + E E h k£a ■ (A - 2) 

i k i k 

If one assumes rigid levels, which is a reasonable approxima- 
tion in many cases, the errors in the energy differences be- 
tween two groups of bands (i.e., across a band gap) can be 
corrected by adding constant energy shifts to the valence and 
conduction band states: 

£ik —> £ik + AEVb for the valence band, and 
£%k — > £%k + AEqb for the conduction band. 

The sum of the band shifts, A£Vb + AEcb, equals the band 
gap error. The expression for the corrected band energy then 
reads 

VB 

i k 
CB 

+ ^2^2Me lk + AE CB ). (A3) 

i k 

Taking the difference between Eqs. JA.2t and (|A.3t , one ob- 



tains the band energy correction term 

AE c b ° rr = (E b - E b ) - (E b - E b ) 

\ / defect \ / ideal 

VB CB 

= ae yb J2 E - x ) +AEcs EE/*' 

i k i k 

Azh Az e 

(A.4) 

where Az/, is simply the number of unoccupied states in the 
valence band and Az e is the number of occupied states in the 
conduction band. 

According to Eq. (01, the defect formation energy of 
charged defects further depends on the position of the valence 
band maximum. The total correction term for the formation 
energy thus reads 

AE corr = , q + A ^) A£ ; VB + Az e AEcB- © 

The energy offsets A£Vb and AEqb can be obtained, for 
example, from GW calculations, which for many systems pro- 
vide band gaps and structures in good agreement with experi- 
ment. 

It is important to realize that this scheme neglects both level 
relaxations and changes in the double counting term. If these 
limitations are acceptable, this method offers a simple a pos- 
teriori correction of the formation energies. It should fur- 
thermore be noted that transition levels are not affected by 
the relative weight of AEVb and A_Ecb- Upon application of 
this correction scheme, acceptor levels track the valence band 
maximum while donor transitions follow the conduction band 
minimum. 
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